Field Induced Reduction of the Low Temperature Superfiuid Density in YBa 2 Cu 3 06.95 
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A novel high magnetic field (8 T) spectrometer for muon 
spin rotation (pSR) has been used to measure the tempera- 
ture dependence of the in-plane magnetic penetration depth 
\ a b in YBa2Cu306.95- At low H and low T, X a b exhibits the 
characteristic linear T-dependence associated with the energy 
gap of a d x 2_ y 2-w&ve superconductor. However, at higher 
fields Xab is essentially temperature independent at low T. 
We discuss possible interpretations of this surprising new fea- 
ture in the low-energy excitation spectrum. 

74.25.Nf, 74.72.Bk, 76.75.-H 



In a superconductor, the resistance to the flow of elec- 
tric current drops to an unmeasurably small value below 
a certain critical temperature T c . This remarkable char- 
acteristic is due to the formation of pairs of electrons (or 
holes), called "Cooper pairs", which link together and 
carry the charge through the sample with virtually no 
opposition. To break apart the pairs, an additional en- 
ergy is needed to excite individual electrons above an 
energy gap which exists at the Fermi surface in the su- 
perconducting state. The nature of these elementary ex- 
citations, known as "quasiparticles" (QPs), is directly 
related to the size and symmetry of the energy gap. The 
gap itself reflects the symmetry of the pair wave function 
(or order parameter) , knowledge of which is essential to 
understanding the physics of the underlying mechanism 
responsible for superconductivity. 

A major breakthrough in the study of high-T c cuprate 
superconductors (HTSCs) came when it was realized that 
the symmetry of the energy gap was different from that 
in conventional low-T c materials. In particular, the en- 
ergy gap was found to vanish along certain directions 
in momentum space. These so-called "nodes" serve as 
a conduit for extreme low-energy QP excitations. One 
of the key early experiments providing evidence for the 
existence of gap nodes was microwave measurements by 
Hardy et al. [fij of the in-plane penetration depth change 
AA a 6 = A a b(T)— A a (,(1.35 K) in the Meissner state of high- 
purity YBa2Cu306.95- In this phase, magnetic field is 
partially screened from the interior by "supercurrents" 
circulating around the sample perimeter. These super- 
currents constitute the response of the superconductor 
to the applied field. The penetration depth A is the char- 
acteristic length scale over which the field decays in from 



the surface, and the quantity A~ 2 is proportional to the 
density of Cooper pairs, i.e. "superfiuid density", n s . 
Because thermal energy can excite QPs, A~ 2 decreases 
with increasing T. In a conventional superconductor, this 
temperature dependence is typically weak at low T be- 
cause the isotropic energy gap exponentially cuts off the 
QP excitations as T^O K. In Ref. M, however, A~ 2 was 
found to decrease sharply upon raising the temperature 
above 1.35 K — the lowest temperature reached in the ex- 
periment. This suggested that the minimum gap size was 
very small. Moreover, at low temperatures AA a {, was ob- 
served to be proportional to T, which is characteristic of a 
superconducting order parameter which has d x 2 _ y 2-w&ve 
symmetry (rather than s-wave symmetry as in conven- 
tional superconductors). Later the same behavior for A a b 
was observed in the vortex state of YBa2Cu2C>6.95, using 
the muon spin rotation (/iSR) technique |2| . In the vortex 
state, magnetic field penetrates the sample in the form 
of quantized flux lines, called "vortices", which usually 
arrange themselves into a periodic array. Screening cur- 
rents circulate around the individual flux lines, so that 
here \ a b is associated with the decay of the field outside 
the vortex cores. 

Recently, researchers have been wondering whether the 
story of the gap symmetry in the high-T c materials is 
complete. Some experiments performed at high mag- 
netic fields in the vortex state show features which are 
seemingly inconsistent with a pure d x 2_ y 2-wetve order pa- 
rameter. For instance, scanning tunneling spectroscopy 
(STS) measurements in YBa2Cu3C>7_,5 (YBCO) suggest 
the existence of localized QP states in the vortex cores [f| , 
which can only exist if there is an energy gap over the en- 
tire Fermi surface. Equally intriguing have been reports 
of an unusual plateau-like feature in the field dependence 
of the thermal conductivity n{H) in Bi2Sr2CaCu20s 
(BSCCO) §-§], and underdoped YBa 2 Cu 3 6 . 63 §. A 
field-induced transition from a d x 2_ J/ 2-wave state to a 
fully gapped state, such as d x 2_ y 2 + id xy or d x 2_ y 2 + is, 
has been offered as a possible explanation for these re- 
sults [pU|]. However, subsequent measurements of k(H) 
have established that the plateau is hysteretic j| and is 
not observed in all samples HJ7]] . These observations sug- 
gest that the vortex lattice (VL) and/or impurities are 
partially responsible for this feature. 

Recent developments in spectrometer design at the 
TRI-University Meson Facility (TRIUMF) in Vancou- 



ver, Canada have made it possible to extend the appli- 
cation of the /iSR technique to magnetic fields as high 
as 8 T. This has allowed us to measure A at fields where 
the anomalies in the STS and n(H) measurements have 
been observed. In a fj.SH experiment the inhomogeneous 
magnetic field associated with the VL is sensitively mea- 
sured by implanting muons into the sample, where their 
spins precess with a frequency which is directly propor- 
tional to the local field. The muon-spin precession signal 
which is obtained by detecting a muon's positron-decay 
pattern, contains the precession frequencies from all the 
muon stopping sites. The line width of the corresponding 
field distribution is roughly proportional to A~ 2 . 

Our measurements were performed on detwinned single 
crystals of YBa2Cu306. 95 (T c = 93.2 K), where a linear In- 
dependence has been previously observed at low H pLllOf . 
To generate a well ordered VL, the sample was cooled in 
the applied field. Neutron scattering measurements on a 
detwinned crystal of YBCO over the field range if = 0.2- 
4 T, show that at low T the vortices arrange themselves in 
a hexagonal lattice, distorted by a-b anistropy |ll[ . This 
is consistent with recent theoretical calculations which 
show that in a d 2 .2„ J/ 2-wave superconductor the VL is 
hexagonal at reduced fields H/H C 2<0.15 Q. 

An example of the measured muon spin precession sig- 
nal is shown in Fig. 1(a). Although the time spectrum 
itself is not very revealing to the eye, a fast Fourier trans- 
form (FFT) can be performed, as in Fig. 1(b), to illus- 
trate the distribution of precession frequencies (or local 
fields) in the sample. We stress that the frequency spec- 
trum in Fig. 1(b) is only an approximation of the ac- 
tual internal field distribution, because the finite time 
range and the "apodization" needed to eliminate "ring- 
ing" in the FFT artificially broadens the output spec- 
trum iQ. We have previously shown that a sharp field 
distribution exists in high quality YBa2Cu30g.95 crys- 
tals, with the basic features slightly smeared out by a 
small random displacement of the vortex positions and 
a small nuclear dipolar contribution p0| . Despite the 
shortcomings of the FFT procedure, the basic features 
of a hexagonal VL are recognizable in the observed aym- 
metric line shape — namely, the minimum field at the cen- 
ter of the triangle formed by three adjacent vortices is 
close to the cusp associated with the Van Hove singu- 
larity produced by the field at the center position be- 
tween nearest-neighbor vortices. The observation of a 
high-field "tail", related to the field in the vicinity of 
the vortex cores, is a signature of a well-ordered VL. 
The inset of Fig. 1(b) shows the field dependence of the 
skewness parameter, a — {(AB) 3 ) 1 / 3 / {{AB) 2 ) 1 / 2 [where 
{(AB) n ) = ((B-(B)) n )}, which characterizes the sym- 
metry of the field distribution. The decrease in a with 
increasing H is predicted to occur in both s- and d x 2_ y 2- 
wave superconductors ]l2j , and does not imply a change 
in VL geometry. The qualitative agreement between 
the behavior of a(H) extrapolated to T = K and at 
T = 50 K, suggests that the VL remains hexagonal over 



this temperature range. 

The muon-spin precession signals were fit assuming an 
analytical Ginzburg-Landau (GL) model for the internal 
field profile associated with the VL Ej] 
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where u 2 = 2^ b G 2 (l + 6 4 )[1 - 26(1 - b) 2 }, B is the av- 
erage internal field, G are the reciprocal lattice vectors, 
^ab is the in-plane GL coherence length and K\{u) is a 
modified Bessel function. The fitting procedure is de- 
scribed in detail elsewhere [Q. We have verified that 
the qualitative results described in this Letter are robust 
with respect to the theoretical model assumed for B(r). 
The solid curves in Fig. 1(a) are an example of a fit of 
the time spectrum to a theoretical complex muon polar- 
ization function P(t) — P x (t)+P y (t), assuming the field 
profile in Eq. ([j]) and a Gaussian distribution of fields for 
the small background signal barely visible in the FFT at 
813.4 MHz. 

Figure 2 shows the temperature dependence of X~ b at 
H = 0.5, 4 and 6 T. The solid curve in Fig. 2, which agrees 
very well with the H = 0.5 T data at low T, represents 
the microwave measurements of AA a & at zero DC field 
|lj , converted to X~ 2 using the fj.SK value of X a b at T = 
1.35 K. At the higher fields the linear T-dependence at 
low temperatures vanishes and the magnitude of X~ b oc n s 
is reduced. This occurs below T«25 K and 35 K at H = 
4 T and 6 T, respectively. While the nature of the QP 
excitations at low T is strongly influenced by the applied 
field, at intermediate temperatures (35 K<T<60 K) the 
field has a negligible effect. At higher T the VL melts 
at a field-dependent temperature T m (H), as shown by 
magnetization measurements on similar crystals [[L6|. In 
the melted region it is not straightforward to extract \ a b- 
The observed decrease in X~ 2 above T m (H) arises from 
a loss of asymmetry in the measured field distribution. 
We now proceed to discuss the possible sources of the 
field-induced reduction of X~ 2 at low T. 

The application of a magnetic field results in a shift 
of the QP energy levels due to the superflow Q . In a 
<i x 2_ a 2-wave superconductor the low-lying QP levels near 
the nodes are shifted below the Fermi surface, resulting in 
a QP current which flows in opposite direction to the su- 
perflow. This leads to a nonlinear relationship between 
the supercurrent density J s and the superfluid velocity 
v s . The net result is a weakened supercurrent response, 
leading to an increased penetration of the field. In the 
Meissner state, Yip and Sauls |lg] predicted that non- 
linear effects in a d x 2_ y 2-w&ve superconductor give rise 
to a linear iJ-dependence for A at T = K. Although 
recent measurements of X(H) Jh],|2(]] show evidence for 
some sort of nonlinear Meissner effect, they do not agree 
with the Yip and Sauls prediction, and do not show a 
saturation at low T as observed here. 

According to Amin, Affleck and Franz [ pl| , a more 
dominant contribution to the field dependence of X a b 



measured by /uSR, comes from the nonlocality of the su- 
percurrent response in the vicinity of the gap nodes. In a 
clean d x 2_ y 2-wa,ve superconductor the coherence length 
£o, which is inversely proportional to the gap size, di- 
verges at the nodes. Thus, near these regions of the 
Fermi surface, £o > A, and J s at a point r is obtained by 
averaging the field over a surrounding region of radius 
£o- For a non-uniform field, nonlocal effects weaken the 



supercurrent response. In Ref. |21|, the effective X a b(H) 
measured by /iSR was calculated for the combination of 
nonlocal and nonlinear effects in a pure d x 2_ y 2-wa,ve su- 
perconductor in the vortex state at T = K. The result 
is that X a b(H) is a nonlinear function of magnetic field, 
predominantly due to a modification of the field distri- 
bution by nonlocal effects. In Fig. 3 we show that our 
measurements of X a b(H) agree extremely well with this 
prediction. Although no prediction for the temperature 
dependence of X a b appears in Ref. pjj , recent calcula- 
tions show a continuous transition as a function of T to 
the T-independent region 122] . We remark that nonlocal 
effects are not expected to play a role in the high-field 
anomalies observed in the STS and thermal conduc- 
tivity pi experiments. 

Another possible explanation for the saturation of X~ b 
at low T is that a field-induced transition to a fully 
gapped state has taken place which freezes out the QP 
excitations. This does not necessarily lead to a change 
in VL geometry — since as noted earlier, the VL is pre- 
dicted to be hexagonal at low H/H C 2 irrespective of the 
pairing symmetry |12|. Furthermore, at low T the vor- 
tices in YBCO are strongly pinned and probably reluc- 
tant to redistribute themselves. However, if a small imag- 
inary component (s- or d xy -wetve) is induced by the field 
at low T, naively one would expect X~ b (T) to saturate 
at a value equivalent to X~ b 2 (T -> K) at H = 0.5 T, 
which is clearly not what we observe. Recently, Yasui 
and Kita 123] have shown by self-consistently solving the 
Bogoliubov-de-Gennes (BdG) equations, that the mixing 
of an s- or d xy -w&ve component, if it exists, must grad- 
ually develop as a function of field in the vortex state. 
Although this is not inconsistent with our /iSR mea- 
surements, the high-field plateau observed in the k(H) 
measurements B indicates a sharp field-induced tran- 
sition of unknown origin. In Ref. [E3J it is also shown 
that the double-peak structure observed near zero bias in 
the STS experiment of Maggio-Aprile et al. || can orig- 
inate from low-energy QP hopping between the vortex 
cores of a pure d x 2_ y 2-w&ve superconductor — negating 
the need for the large d xy component suggested by Franz 
and Tesanovic B in an earlier treatment of an isolated 
vortex using the BdG formalism. Thus, current theories 
can explain both the STS and /iSR anomalies at high field 
in terms of a pure d x 2_ y 2-w&ve superconductor. We note 
that while the peak structure was not observed in STS 
measurements of the vortex cores in BSCCO by Renner 
et al. |24] , a recent STS study by Davis et al. p5[ shows 
evidence for QP core states in this material. 



Finally, Fig. 4 shows the field dependence of the in- 
plane coherence length £ a b from the fits to the /iSK spec- 
tra, extrapolated to T = K. Previously, we established 
that £ a {, is nearly proportional to the size of the vortex 
cores ro, and that tq expands at low fields [E6J. The 
high-field measurements clearly show that the core size 
saturates with £ t>~ 18.5 A when the density of vortices 
(~ H) becomes large. The expansion of ro at low fields 
is qualitatively consistent with the predictions from a re- 
cent quasiclassical theoretical study of the vortex struc- 
ture in pure s- and d x 2_ y 2-w&ve superconductors |12j. 

In summary, we observe a field-induced saturation in 
the temperature dependence of X a b at low T measured by 
/iSR. The behavior of X~ b (T) appears inconsistent with 
a transition to a mixed d x 2_ y 2+id xy (d x 2_ y 2 +is) state. 
The measured field dependence of X a t is in agreement 
with the theory of nonlocal and nonlinear effects in the 
vortex state of a pure d x 2_ y 2-wave superconductor. 
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FIGURE CAPTIONS [18] 

Figure 1. (a) The muon spin precession signal at [19] 
T = 20 K and H = 6.0 T displayed in a reference frame 

rotating at about 3 MHz below the Larmor precession [20] 

frequency of a free muon in the external field [Note: A [21] 
is the maximum precession amplitude], (b) The FFT of 

(a) using a Gaussian apodization with a 3 /xs time con- [22] 

stant. Inset: field dependence of the skewness parameter [23] 

a extrapolated to T = K [solid circles] and at T = 50 K I 24 ] 

]open circles]. !■ J 

L J [26] 

Figure 2. Temperature dependence of X~ b 2 at H — 0.5 T 
, 4 T and 6 T. The solid curve represents the zero-field 
microwave measurements of AA a (,(T) in Ref. M. 

Figure 3. Magnetic field dependence of X a b hi 
YBa 2 Cu306.95 extrapolated to T — K [open circles] 
and the predicted behavior at T — K from Ref. |EM for 
the combination of nonlinear and nonlocal effects in the 
vortex state of a d x 2 _ y 2-w&ve superconductor [solid dia- 
monds] . 
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